Interference oscillations of microwave photoresistance in double quantum wells 
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We observe oscillatory magnetoresistance in double quantum wells under microwave irradiation. 
The results are explained in terms of the influence of subband coupling on the frequency-dependent 
photoinduced part of the electron distribution function. As a consequence, the magnetoresistance 
demonstrates the interference of magneto-intersubband oscillations and conventional microwave- 
induced resistance oscillations. 
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The phenomenon of microwave-induced resistance 
oscillations 1-3 (MIRO) in two-dimensional (2D) elec- 
tron systems under perpendicular magnetic fields has at- 
tracted much interest. 4 These oscillations are periodic 
in the inverse magnetic field with a period determined 
by the ratio of the microwave radiation frequency lu to 
the cyclotron frequency u c and survive at high temper- 
atures. Basically, the observed oscillatory photoconduc- 
tivity is caused by the Landau quantization of electron 
states, though different microscopic mechanisms of this 
phenomenon are still under discussion. 

The influence of microwave irradiation on the magne- 
totransport properties is currently under investigation in 
quantum wells with a single occupied subband. In our pa- 
per, we underline the importance of similar studies for the 
systems with two occupied 2D subbands, where the mag- 
netotransport shows special new features as compared 
to the single-subband case. Apart from the commonly 
known Shubnikov-de Haas oscillations (SdHO), there ex- 
ist the magneto-intersubband (MIS) oscillations of resis- 
tivity caused by periodic modulation of the probability 
of intersubband transitions by the magnetic field 5,6 (see 
Ref. 7 for more references) . These oscillations survive at 
high temperatures, because they are not related to the 
position of the Landau levels with respect to the Fermi 
surface. 

Our recent observation 7 of large-amplitude MIS oscilla- 
tions at magnetic fields below 1 T in high-mobility double 
quantum wells (DQWs) has established that these oscilla- 
tions are a well-reproducible feature of magnetotransport 
in such systems. The present measurements, supported 
by a theoretical analysis, suggest that the behavior of 
oscillating magnetoresistance in DQWs under microwave 
irradiation is caused by an interference of the physical 
mechanisms responsible for the MIS oscillations and con- 



ventional MIRO. We have studied dependence of the re- 
sistance of symmetric balanced GaAs DQWs with wells 
widths of 14 nm and different barrier widths df, = 1.4, 
2, and 3 nm on the magnetic field B in the presence of 
microwave irradiation of different frequencies and at dif- 
ferent temperatures and intensities of radiation. 
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FIG. 1: (Color online) Measured (upper panel) and calcu- 
lated (lower panel) magnetoresistance of DQWs for different 
intensities of 170 GHz radiation at T = 1.4 K. The effect of 
the radiation can be neglected for the attenuation of 60 dB. 



2 



The samples have high mobility of 10 6 cm 2 /V s and 
high total sheet electron density n s ~ 10 12 cm -2 . The 
samples were mounted in a waveguide with different cross 
sections. The measurements were performed for perpen- 
dicular and parallel orientation of the current with re- 
spect to microwave polarization. The resistance R = R xx 
was measured by using the standard low-frequency lock- 
in technique. No polarization dependence of the resis- 
tance has been observed. While similar results were ob- 
tained for samples with different barrier widths, the data 
reported here corresponds to 1.4 nm barrier samples. 
Without the irradiation, the magnetoresistance shows 
large-period MIS oscillations clearly visible starting from 
B = 0.1 T, and it also shows small period SdHO super- 
imposed on the MIS oscillation pattern at higher fields 7 
(see Fig. 1). The subband separation in our structure, 
A12 = 3.67 meV, has been found from the MIS oscillation 
periodicity at low B. 
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FIG. 2: (Color online) Measured (upper panel) and calcu- 
lated (lower panel) magnetoresistance of DQWs under 170 
GHz irradiation for different temperatures. The inset shows 
the dependence of electron temperature T e on the magnetic 
field. Weak oscillations of T e are caused by oscillations of the 
absorbed radiation power. 

An increase in the intensity of the radiation leads to 
the expected damping of the SdHO owing to electron 
heating, and also causes suppression of some groups of 



MIS peaks, followed by the inversion (flip) of these peaks 
at high intensity, while the other MIS peaks are con- 
tinuously enhanced by the irradiation. With increasing 
temperature, the inverted peaks go up and the enhanced 
peaks go down, so the conventional MIS oscillation pic- 
ture is restored (Fig. 2). Measurements at different fre- 
quencies (Fig. 3) have confirmed that the selective peak 
flip is correlated with the radiation frequency and follows 
the periodicity determined by the ratio uj/lu c . This al- 
lows us to attribute the observed effect to MIRO-related 
phenomena. However, the peak flip cannot be explained 
by a simple superposition of the factors cos(27rAi2//iw c ) 
and — sin(27ra>/w c ) describing the MIS oscillations and 
the MIRO, respectively. 
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FIG. 3: (Color online) Measured (upper panel) and calculated 
(lower panel) magnetoresistance of DQWs at T — 1.4 K for 
different frequencies of microwave excitation. The curves for 
100 and 170 GHz are shifted up for clarity. 

The theoretical explanation of our data is based on 
the physical model of Dmitriev et al., s generalized to 
the two-subband case and improved by taking into ac- 
count electrodynamic effects in the absorption of radia- 
tion by 2D layers. 9-11 We assume that the elastic scat- 
tering of electrons, including the intersubband scatter- 
ing, is much stronger than the inelastic one, so the elec- 
tron distribution function (averaged over the time pe- 
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riod 2tt/lo) under microwave excitation remains quasi- 
isotropic and common for both subbands, though essen- 
tially non-equilibrium. This energy distribution function, 
f e , is found from the kinetic equation G £ (f) = —J e {f), 
where G e is the generation term, and J e is the inelastic 
collision integral. Under the excitation by linearly polar- 
ized high-frequency electric field E t = Ecos(wi), weak 
enough to neglect the multi-photon processes, the gener- 
ation term is written as 



G F 



1 



D F 8irmuj 2 



+d>)(/e-^-/e)j. (1) 

where D £ is the density of states, e is the electron 
charge, and m is the effective mass of electron. Next, 
= J BA£ ±) , where a£ ±) = [1 + 2na ± {uj) / c^fe]- 1 , are 
the amplitudes of the circularly polarized components of 
electric field in the 2D plane, found from the Maxwell 
equations. 9,10 Here, cr±(cu) are the complex conductiv- 
ities describing response of the electron system to the 
circularly polarized fields, c is the velocity of light, and e 
is the dielectric permittivity of the medium surrounding 
the quantum well. The functions <f>E^ (lu), which describe 
the probability of electron transitions between the states 
with energies e and e+fku, are written in terms of Green's 
functions as ^^(lu) = Re(Qf^ ± " > — Q£,v^), where 

Qt i+) = %T E E >/(" + D(»' + 1) E 
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x ((G» > s (n + l Py ,n' + lp' v )Gi^(n'P>Pv))) > (2) 

and QefJ" is given by the permutation of the indices, 
n + 1 <-> n and n' + 1 <-» n' , in the arguments of Green's 
functions in Eq. (2). The Green's functions G, retarded 
(R) and advanced (A), are determined by interaction of 
electrons with static disorder potential in the presence of 
magnetic field. They are written in the representation 
given by the product of the ket-vectors \j) and \np y ), 
describing the 2D subbands and the Landau eigenstates, 
respectively. Here j = 1,2 numbers the electron subband 
of the quantum well, n is the Landau level number, and 
p y is the continuous momentum (the Landau gauge is 
used). The double angular brackets in Eq. (2) denote 
random potential averaging. The conductivities <j±{lo) 
are expressed in terms of the functions (2) as 
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In the limit w — ► 0, Eq. (3) gives the dc conductivity 
components a d — Re(cr + + cr_)/2 and = Im(cr + — 
a-)/2. 



Analytical evaluation of the correlation functions in 
Eq. (2) can be done in the case of relatively weak mag- 
netic fields, by using the self-consistent Born approxima- 
tion and by expanding the Green's functions in powers of 
small Dingle factors. In application to balanced DQWs, 
where the subband energy separation is typically small 
compared to the Fermi energy, we assume that the dif- 
ference in subband populations is small compared to n s . 
It is justified for our samples where the Fermi energy 
(17 meV) is much larger than A 12 /2 = 1.8 meV. The 
difference between quantum lifetimes of electrons in the 
two subbands, t\ and T2, (and between the correspond- 
ing transport times) is also small, so we neglect it ev- 
erywhere except the Dingle exponents. In the first order 
in the Dingle factors dj = cxp(— 7r/w c Tj), the generation 
term acquires the following form: 
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s± = (ui±Lj c ) 2 T 2 r7 and T tr is the transport time, common 
for both subbands in our approximation. 

The kinetic equation is solved by representing the dis- 
tribution function as a sum /° + Sf e , where the first 
term slowly varies on the scale of cyclotron energy, while 
the second one rapidly oscillates with energy. 8 The first 
term satisfies the equation ^r£ 0) (f° +hw + /°_ ftw - 2/ e ) = 
— J e (/°). Assuming that the electron-electron scatter- 
ing controls the electron distribution, one can approx- 
imate /° by a heated Fermi distribution, /° = {1 + 
cxp[(e — £ir)/T e ]} _1 . The electron temperature T e is 
found from the balance equation P u = P p h 1 where P u = 



JdeeD E G E (f) = Re[|AL + T^+M + |AiT'|V-M]E 2 /4 
is the power absorbed by the electron system and P p h = 
— J ds sD £ J e (f°) is the power lost to phonons. The sec- 
ond term satisfies the equation 



r^[5f E+hw +Sf e . hw -25f s ] 
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where we have used the relaxation time approximation 
for the collision integral, J e (Sf) ~ —6f £ /Ti n , with inelas- 
tic relaxation time Tj„. Solution of this equation is 



Sfe ~ — 7^— > cLsin — ^- — , (8) 
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where P u = AQ^rl >T in . 

Using the energy distribution function found above, 
one can calculate the dc resistivity p xx = 0^/(0^ + 0jj: 

^ = 1 - 2%g £F + ^(d 2 + d\) + 2^d 1 d 2 

2ttAi2 A u ( l9 , , 27rAi 2 \ , , 

x cos \d\ + d\ + 2dxd 2 cos - , (10) 

nu! c 2 \ hLo c ) 

where po = m/e 2 n s T tr is the zero-field Drude resistiv- 
ity and A12 = £2 — £1 is the subband separation. The 
second term, proportional to g SF , describes the SdHO. 
The third and the fourth terms describe positive magnc- 
toresistance and the MIS oscillations, 7 respectively; these 
terms are written under condition of classically strong 
magnetic fields, u c T tr S> 1. Here, t{\ = r|5 and t% are 
the intrasubband and intersubband transport scattering 
times, which contribute to the total transport time ac- 
cording to l/r tr = 1/t*i + 1/ T i2- The term proportional 
to the oscillating factor describes modification of the 
oscillatory resistivity under photoexcitation. Another ef- 
fect of the excitation is the electron heating, which in- 
creases thermal suppression of the SdHO, described by 
the factor T e = (2i: 2 T e /huj c )/ smh(2n 2 T e /fuv c ). 

The most essential feature of the resistivity given by 
Eq. (10) is the presence of the product of the oscil- 
lating factors cos(27tAi2/?icl | c ) and sin(27rw/a; c ), which 
corresponds to an interference of the MIS oscillations 
with the MIRO. The interference takes place because 
the photoexcitation involves both subbands [see Eqs. (8) 
and (9)], so the scattering- induced coupling between the 
subbands, which oscillates as cos(27rAi 2 /fiu c ), modifies 
the amplitude of the oscillating photocurrent. Usually, 
the microwave quantum energy is smaller than the sub- 
band separation, and the magnetoresistance shows fast 
MIS oscillations modulated by a slow MIRO component 
oc — sin(27rw/oj c ). 

In our calculations, we have used Eq. (10) with 
di = d<i = exp(— 7t/cj c t) and t\ 2 = t{\ = 2r tr , which 
is a good approximation for balanced DQWs. 7 We have 
taken into account heating of electrons by the field and 
dependence of the characteristic times r, Tj r , and Tj„ 
on the electron temperature T e . The latter was deter- 
mined by using the collision integral for interaction of 
electrons with both deformation and piezoelectric poten- 
tials of acoustic phonons. The inelastic scattering time 
Ti n is assumed to scale with the temperature as T~ 2 , 
with H/Tin = 4 mK at T e = IK, according to the 
theoretical estimates for electron-electron scattering 8 ap- 
plied to our samples. The temperature dependence of the 
transport time r tr and quantum lifetime r of electrons in 



our samples has been empirically determined by analyz- 
ing temperature dependence of zero-field resistance and 
MIS oscillation amplitudes. The microwave electric field 
E ~ 2.5 V/cm is estimated by fitting calculated ampli- 
tudes of magnetoresistance oscillations to experimental 
data; the corresponding electron heating is in agreement 
with the observed suppression of the SdHO amplitudes. 

The results of the calculations, shown in the lower 
parts of Figs. 1-3, are in good agreement with ex- 
periment, and capture all qualitative features of the re- 
sistance dependence on the magnetic field, temperature, 
radiation power, and frequency. Note that the heating 
of electrons by microwaves is not strong at lo c ~ u> (sec 
inset to Fig. 2) because of radiative broadening of the cy- 
clotron resonance 10 ' 11 in our samples. The only feature 
that is not described by the theory is the observed re- 
duction in the MIS oscillation frequency 7 at B > 0.25 T, 
which is possibly related to a decrease in the subband sep- 
aration owing to enhanced Coulomb correlations and/or 
modification of screening in magnetic fields. 

Understanding the MIRO physics requires comparison 
of measured photoresistance with the results provided by 
existing theories. Among the physical mechanisms 12 re- 
sponsible for the MIRO, the inelastic mechanism, 8 ' 12 de- 
scribing the oscillations as a result of microwave-induced 
change of the isotropic part of electron distribution func- 
tion, is predicted to dominate at moderate radiation 
power, owing to a large ratio of T in /r, which is of the 
order 10 2 under typical experimental conditions, includ- 
ing our experiment. The observed insensitivity of the 
MIRO to the polarization of incident radiation, and T~ 2 
scaling of the oscillatory photoresistance amplitude 11 are 
in favor of the inelastic mechanism. We have shown that 
application of the basic principles of the inelastic mech- 
anism theory 8 to the systems with two-subband occupa- 
tion confirms the reliability of theoretical estimates for 
the inelastic relaxation time and leads to a satisfactory 
explanation of the new features reported in this paper. 

In conclusion, we have presented experimental and 
theoretical studies of oscillatory magnetoresistance in 
DQWs under microwave irradiation. The new kind of 
interference oscillations observed in these systems ap- 
pears because the photoinduced part of the electron dis- 
tribution, which oscillates as a function of microwave fre- 
quency, is modified owing to subband coupling and be- 
comes also an oscillating function of the subband sepa- 
ration. 
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